Quantum Tunneling Detection of Two-photon and Two-electron Processes 
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We analyze the operation of a quantum tunneling detector coupled to a coherent conductor. We 
demonstrate that in a certain energy range the output of the detector is determined by two-photon 
processes, two-electron processes and the interference of the two. We show how the individual 
contributions of these processes can be resolved in experiments. 



The quantum nature of electron transfer in coherent 
conductors is seldom explicitly manifested in averaged 
current-voltage curves. To reveal it one should measure 
current noise and/or higher-order correlations of current 
comprising Full Counting Statistics which arise from the 
transfer [1| . Such measurements not only reveal the dis- 
crete nature of the charges transferred, they also quan- 
tify quantum many-body effects in electron transport and 
may be used for the detection of pairwise entanglement of 
transferred particles 0, EJ 0, 0] • If the noise is measured 
at frequencies in the quantum range, Tuo fceT, the 
measurement amounts to the detection of photons pro- 
duced by the current fluctuations. This aspect is impor- 
tant in view of attempts to transfer quantum information 
from electrons to photons and back Q. 

It was demonstrated in 

7] 

that one needs a quantum 
detector to measure quantum noise. Indeed, any classical 
measurement of a fluctuating quantity would give a noise 
spectrum symmetric in frequency, S(u>) = S(—u>). A 
quantum tunneling detector is generally a quantum two- 
level system with a level separation e > 0. The result 
of detection are two transition rates: r up from the lower 
to the higher level and Tdown for the reverse direction. 
The most probable transitions are accompanied by either 
absorption or emission of a photon of matching energy 
Hlu = e. One can define the noise spectrum in such a way 
that it is proportional to the transition rates S(±e/h) oc 
rup,down(£)- Differences between r uPj down thus manifest 
the quantum nature of noise. If the source of noise is 
a coherent conductor biased by a voltage V, detector 
signals in the range e < eV are readily interpreted in 
terms of single electron transfers through the conductor. 
The maximum energy gain available for electrons in the 
course of such transfer is eV. Consequently this value 
also limits the energy of the emitted photon. 

A first proposal for the experimental realization of a 
quantum tunneling detector included transitions between 
two localized electron states in a double quantum dot || . 
However it does not matter much if the tunneling oc- 
curs between localized or delocalized electron states and 
if all tunnel events are accompanied by the same energy 
transfer e. In most practical cases the energy dependence 
of the rates r up < jown can be readily extracted from the 
measurement results. This is why quantum tunneling 



detection has been experimentally realized in a super- 
conducting tunnel junction and in a single quantum 

dot H3. 

In this Letter we study quantum tunneling detection 
in the range eV < e < 2eV assuming e, eV » k^,T. 
The motivation is that for these energies the detec- 
tor is not sensitive to single-electron one-photon pro- 
cesses described above and its output — the transition 
rate r up — is determined by much more interesting two- 
particle processes. It is clear from plain energy consid- 
erations that transitions may originate from two-photon 
processes. Such two-photon absorption can occur given 
any non-equilibrium photon distribution bounded by eV, 
not necessarily produced by a coherent conductor. Less 
obvious and specific for a coherent conductor is a cooper- 
ative two-electron process. Indeed, if two electrons team 
up in crossing the conductor they can emit a single pho- 
ton with an energy up to 2eV. Essential for this cooper- 
ation are electron-electron interactions. It is known [Tl| 
that the most important electron-electron interaction in 
this energy range is due to the electromagnetic environ- 
ment of the conductor, the same environment in which 
the non-equilibrium photons dwell. We quantify the sig- 
nals due to two-photon and two-electron events and find 
them to be of the same order of magnitude. We also show 
that part of the signal is due to quantum interference 
of these two processes: one-and-half-photon absorption 
events. We demonstrate how different contributions can 
be separated in experiments thereby facilitating the di- 
rect observation of two-particle processes in the context 
of quantum transport. 







<9(t) 1 








) 


7 


A 




T 

n 




C 



V 



FIG. 1: Model circuit 
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We concentrate on a model circuit consisting of four 
elements as given in figure ^ A voltage biased coherent 
contact characterized by a set of transmission eigenval- 
ues {T n } is embedded in an electromagnetic environment 
with impedance Z u . The environment transforms the 
current fluctuations in the conductor to voltage fluctua- 
tions in node A which are conveniently expressed in terms 
of a phase <p = f J dtV(t). The most general model in- 
cluding the detector and coherent contact would be a 
four-pole circuit studied in [12j that couples two poles of 
the detector with two poles of the contact. Here, we con- 
centrate on the experimentally relevant case of capacitive 
coupling. Owing to voltage division between two capac- 
itors, the detector senses a fraction < a < 1 of the 
voltage fluctuations in node A. We will see that chang- 
ing the "visibility" parameter a, enables the separation 
of two-electron and two-photon processes in experiments. 
The relevant impedance is made up of an environmental 
impedance combined with that of two capacitors and that 
of the coherent contact. We measure this impedance z u 
in units of Rk = 2irh/e 2 and assume the low-impedance 
limit z u « 1; this provides us with a physically justified 
small parameter. 

The detector consists of two localized charge states 
connected by a tunnel amplitude T. In the presence of 
voltage fluctuations in the node A, the amplitude is mod- 
ified as follows: T(t) -^Te lav ^ t '. In perturbation theory, 
the inelastic tunneling rate between two states separated 
by e is given by correlators of oup{t) 0] 



2irk 2 



dt(e 



ioup(t) e -ioup{0)\ e ^et 



(1) 



The rate T(e) is therefore the Fourier transform of the 
correlation function (e^We-^ )}, aip(t) being the 
phase fluctuations over the detector. From now on we 
take h = e = fcs = 1- 

Equation tells us that the inelastic tunneling rates 
in the detector are completely determined by the voltage 
fluctuations over the junction. Therefore measuring the 
inelastic current through the dots we are sensitive to the 
noise spectrum of the environment. 

To evaluate ( e <Q! ¥>(*) e -"*v(o)) we construct a path 
integral representation of this quantity using a non- 
equilibrium Kcldysh technique for quantum- 
circuits [iH 



a iaip(t) —ia<p(0) 



T>[4>] exp{-iS eBV [<p] 



iS cond [cf)} + ia[-(p + (0) + <p~(t)]}. 



(2) 



The integration goes over the time-dependent fluctuating 
fields ^5 ± (i) in node A, ± corresponding to the forward 
(backward) part of the Keldysh contour. S cnv and S con d 
are the contributions to the Keldysh action originating 
from the environment and the coherent conductor respec- 
tively. 



Since the environment is linear, its action is quadratic 
in the fields and at zero temperature reads (cf. |l5j) 



Scnv = J duJ <f> T LJ A(ld)<t> u 



(3) 



with 



Zu, being the corresponding impedance. We use Fourier 
transformed fields <p u = (0 w ,Xw) T defined with ip ± = 

All non-quadratic contributions to the action originate 
from the coherent conductor. The action S COVl & can be 
expressed in terms of Keldysh Green functions Gl,r of 
electrons in the reservoirs left and right of the contact 


S cond = ~ £ Tr ln[l + %{G L (0), Gr} - 2)]. (4) 



The fields <j> enter in this action via the gauge transform 

To comprehend the physics involved, let us first disre- 
gard any non-quadratic parts and take only the quadratic 
part of Scond- In this case the path integral is Gaussian, 
and can be evaluated exactly. We recover the well known 
result from P(£)-theory (cf. HEl): (e^Me-^vC )) = 
exp[J(i)] with 

J(i) = {oap{t)aip$)) =a 2 J dw^-^-K (u)^* - 1]. 

(5) 

The impedance includes the dimensionless conductance 
g c = J2 n Tn of the contact. 

In the limit of T = we find in agreement with results 

of Hi 

K(u) = g c {FD{uj + V) + (2 - 2F)D{u) 

1 (6) 
+ FD(u - V)} + 2Rc{ — }D(uj) 

with D(uj) = —Lo6(—Ld) and the Fano factor F = 
£„T„(1 - T n )/J2 n T n . The first term in K{u) (oc g e ) 
represents the non-equilibrium current noise spectrum of 
the coherent contact that vanishes for u> > V. In physical 
terms this means that the highest energy u> an electron 
can emit traversing the conductor is exactly V. The sec- 
ond part represents the spectrum of the environment. It 
is zero for li > 0, since the environment can only absorb 
energy at T = 0. 

The time-dependent part of J(t) is the Fourier trans- 
form of K{lu)\zuj\ 2 /to 2 and r(e) is in turn the Fourier 
transform of exp[J(t)]. If we expand exp[J(t)] in terms 
of J(t), the n-th term presents the contribution of a 
process involving absorption of n photons in the detec- 
tor. Such an n-photon process dominates in the interval 
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FIG. 2: n-photon contributions (dashed) to the detector out- 
put. Each contribution dominates in the energy interval 
(n — 1)V < e < nV. Subsequent contributions are suppressed 
by a factor « a 2 z 2 g c F. This is seen as a staircase structure in 
the log plot. The sum of all contributions is given by the solid 
line. To produce this plot, we took z^/u — z/V, z = 0.01, 
g c F = 1.75, a = 1. 



l)V < e < nV and its contribution is proportional 



to a 



2n 



The one-photon contribution gives Tup /Fdown — zg c F. 
Each extra photon brings in a small factor, such that 
r("+!)/r(™) — a 2 z 2 g c F. This is seen as a staircase in 
the log plot presented in Figure |3 

What we did was wrong since we did not take into 
account the non-quadratic terms in the action. These 
describe more interesting many-electron processes and 
are — as we show below — of the same order of magnitude. 
Since the path integral in @ can not be evaluated in 
general, we proceed by perturbative expansion. 

Indeed, since \z\ -C 1, the Gaussian part of the action, 
being proportional to z , suppresses fluctuations in the 
path integral and we can treat the remaining part pertur- 
batively. First we expand — iS con d[(f>\ around <f> — 0. As 
mentioned previously (see the discussion below equation 
©), the first and second order terms just renormalize 
the impedance. The exponential of the remaining higher 
order terms is then again expanded in <p around = 0. 
This expansion may be represented in terms of diagrams 
such as those in figure |3J Diagram (a) represents a high 
order term, from which the general structure becomes 
clear: Diagrams consist of lines, polygons and external 
vertices. The expansion contains not only connected dia- 
grams, but all disconnected diagrams as well. A polygon 
with n vertices is associated with the symmetrized n-th 
order coefficient in the Taylor expansion of — iS con d[cf>\- 
Each polygon contributes a factor g c . Lines represent 
propagators of <fi corresponding to the Gaussian action 
with renormalized impedance. They are of order z mak- 
ing n-line diagrams z n in leading order. External vertices 
(indicated by dots in the figure) are associated with the 
time-dependent linear term ia[—ip + (0) + <p~ (t)] in equa- 
tion J2J. Thus a diagram with s external vertices gives 



FIG. 3: Typical Feynman diagrams appearing in the expan- 
sion of the path integral representation of the tunneling rate. 



a correction proportional to a s . Furthermore, diagrams 
without external vertices are time-independent and ac- 
cording to equation Q contribute only to elastic tunnel- 
ing processes. Diagrams (b) to (f) represent some of the 
lowest order terms in the expansion. 

We consider transitions in the detector where energies 
between V and 2V are absorbed. In this interval, the 
detector output is given by diagrams (b), (c) and (d), 
which are proportional to a 4 , a 3 and a 2 respectively. 

From the results presented in figure [2] we have learned 
that n-photon processes come with a coefficient a 2n . 
Hence the a 3 contribution is not readily expected: it 
seems to signal a process with one- and-half photons ab- 
sorbed. We disregard diagram (e) which only contributes 
to elastic processes. In the energy interval considered, 
the combined z 3 contribution of the included diagrams 
is zero and we obtain a tunneling rate that goes as g 2 z 4 . 
Since a diagram like (f) has four lines, it could poten- 
tially contribute to the current with the same order in z. 
However, its contribution can only be proportional to g c 
and is disregarded. 

The expansion of Scond up to fourth order terms and 
subsequent evaluation of the diagrams is straightforward 
but requires rather involved and lengthy calculations. 
Fortunately in the interval of interest the three contri- 
butions can be combined in a compact expression 

r v 

T up = 2\T\ 2 g 2 c F 2 / dw{V ~ w)(e - V — oj) 

Js-V 



2 

a z F 



2 e- 



h a — 

lo e 



(7) 



which is the main result of our work. The rate is propor- 
tional to the square of the zero-frequency current noise 
5 co „d(0) = | g c F. 

The part proportional to a 4 (diagram b) represents a 
two-photon process originating from the quadratic part 
of 5 CO nd and was already present in figure We have 
thus shown that there are contributions of the same or- 
der resulting from non-linearities in the conductor. The 
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a 2 term (diagram d) is the result of the two-electron 
and one-photon process expected from general reason- 
ing presented in the introduction. We see that the a 3 
term comes from the cross-term in the modulus square. 
This unambiguously identifies digram (c) as the result of 
quantum interference of the two-electron process and the 
two-photon processes — an interpretation that was not 
obvious from the beginning. 

To understand this interference, we note that the pho- 
ton modes involved are delocalized across the whole cir- 
cuit. A photon in each mode can be absorbed in the 
detector as well as in the environment or the contact. 
An elementary process is such that the final state dif- 
fers from the initial state by two photons absorbed in 
two given modes. The final state can be reached by two 
amplitudes: one with both photons absorbed in the de- 
tector and one with a photon absorbed in the detector 
and a photon absorbed in the environment. While the 
squares of these amplitudes represent the probabilities of 
two-photon and two-electron processes respectively, their 
cross-term gives rise to an interference contribution oc a 3 . 

The simplest concrete model is that of a frequency- 
independent impedance, z u = z at u> ~ V . The integra- 
tion in Eq. yields for the three distinct contributions ( 
s = e/V, 1<£<2) 



z 4 \T\ 2 q 2 .F 2 



a 

2a 3 
a 2 



ln(e - 1) 



(8) 



-fln(e-l)- 



4(2-£) 



All contributions scale as (e — 2V) 3 at the two-photon 
threshold and logarithmically diverge at approaching the 
one-photon threshold, (see Fig. HJ. 

Here, we have quantified the contributions for a very 
specific non-linear quantum noise source: a coherent con- 
ductor. However in the case of any unknown source of 
this kind the a-dependence of the contributions allows 
one to separate and identify them experimentally (right 
pane of Fig. 0J. One would measure the detector output 
changing the coupling to the detector. Formally, three 
measurements at three different a are sufficient to deter- 
mine the relative strength of all three contributions. In 
any case, in the limit of small coupling a — > the de- 
tector output is dominated by two-electron events. Fur- 
ther characterization may be achieved by engineering of 
a frequency-dependent impedance. For instance, setting 
z(ui — e) to kills both interference and two-electron 
contribution. 

To conclude, we have shown that the quantum tun- 
neling detector in the energy interval specified is selec- 
tively sensitive to two-particle processes. The detector 
output is generally determined by three contributions: 
two-photon processes, two-electron processes and the in- 
terference of the two. These three sources can be dis- 
tinguished experimentally by measuring at different cou- 
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FIG. 4: Left: Contributions to the detector output due to 
two- photon (solid), two-electron (dashed) processes and their 
interference (dash-dotted) at a = 0.8 versus detector level 
splitting. The inset presents a zoom at e — > 2V . Right: 
Different dependence on the coupling strength a enables ex- 
perimental identification of these three contributions. The 
contributions and their sum (dotted) are plotted for e = 1.3V. 



plings a to the detector. Our results thus facilitate the 
direct observation of many-particle events in the context 
of quantum transport. This work was supported by the 
Dutch Foundation for Fundamental Research on Matter 
(FOM). 
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